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Abstract. The photon motion in a Michelson interferometer is re-analyzed in both 
geometrical optics and wave optics. The classical paths of the photons in the 
background of gravitational wave are derived from Fermat principle, which is the same 
as the null geodesics in general relativity. The deformed Maxwell equations and the 
wave equations of electric fields in the background of gravitational wave are presented 
in flat-space approximation. Both methods show that the response of an interferometer 
depends on the frequency of a gravitational wave, however it is almost independent of 
the frequency of the mirror’s vibrations. It implies that the vibrating mirror cannot 
mimic a gravitational wave very well. 
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1. Introduction 


The transient gravitational wave (GW) events, GW150914 and GW151226, were 
detected in the first run (O1, from September 12, 2015 to January 19, 2016) of the 
advanced LIGO [1], [2]. The events were reported to be the results of the coalescence 
of two black holes of masses 36Moand 29Mo to a black hole of 62M. and of two black 
holes of masses about 14M. and 8Mo to a black hole of mass about 21Mo, respectively. 
The GW signals are chirps in the frequency ranges between 35 to 150 Hz and between 35 
to about 700 Hz, respectively. The signals of GW150914 and GW151226 are extracted 
from the sea of noises by the same methods (for example, see refs. [3], [4]). Hence, the 
explanation of the observations are widely accepted. 
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A LIGO detector is, basically, a Michelson intermerometer (MI) with its arms 
replaced by Fabry-Pérot (FP) cavities [5]. One of the roles of the FP cavities is to 
increase the interaction time of photons with GWs. On average, a photon travels 140 
round-trips in the FP cavity. In other words, a photon stays in a cavity for about 
280L/c ~ 3.73 ms on average, where L ~ 4 km is the average length of arms. In 3.73 ms, 
a GW with frequency about 268 Hz will propagate through the detector for a wavelength 
and thus a photon cannot ‘feel’ the length variation of the arms [6], [7]. However, it 
is widely accepted that the response of a Michelson-Fabry-Pérot interferometer is as 
about 2N times as the response of a same size Michelson interferometer, where N is the 
average number of the round trips (see, for example, [8],{3]). The previous studies on 
the instruments and calibration [8], [9], [10], [11], [12], [13], [14] show that the response 
of a detector to a GW will not depend on the frequency when the frequency of a GW 
is lower than 1 kHz. Why is there such a kind of inconsistency? 

In fact, almost all previous studies on the response of detector to a GW are 
conducted by use of the swing of mirrors [9], [10], [11], [13], [14]. They seem to be 
supported by the theoretical analysis for the response to a GW [12]. In these analysis, 
the Laplace transformation is used, which is still questionablet. In the present paper, 
we shall re-analysis the response of a detector to a GW, without the help of Laplace 
transformation. We shall focus on the Michelson interferometer in this paper. 

The motion of photons in a detector is govern by general relativity. In other words, 
the photons should move along the null geodesics in the space-time. In fact, there are 
two other ways to describe the motion of photons. The first one is the method of classical 
geometrical optics, in which paths of photons in three dimensional space are determined 
by the Fermat principle. To apply the Fermat principle in vacuum of a curved spacetime, 
the metric is regarded as a position-, direction-, and time-dependent refractive-index 
tensor of a special transparent medium in a flat space. The second one is the method of 


t To our knowledge, the first application of the Laplace transformation in the study of the response of 
an interferometer is in [15] and [8]. In the references, the following equation is used 


L£ cf naat) = = (ate), 


where £ represents the Laplace transformation, p is the Laplace domain variable. However, for 
h(t) = sin(wt), a direct calculation shows 


z(pcos(wL) + wsin(wL) — p), 
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They are not equal to each other. The equation should be revised as 
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classical wave optics, which is based on the wave equation of electric field derived from 
the Maxwell equations. The first aim of the present paper is to study the motion of 
photons in a single arm by use of the two methods. Although a realistic detector needs 
double arms, we may just consider a single one because the motion of photons in each 
arm is govern by the same law but has the different polarity due to the quadruple nature 
of the GW. We shall show that the paths of photons determined by the Fermat principle 
are the same as that determined by the null geodesics even in the background of GW. 
We shall also show that a GW may serve as a time-dependent medium in flat space-time 
and present the deformed Maxwell equations and wave equations of electric field in the 
flat space-time approximation. With the two methods, we may obtain the same claim 
as that by use of the null geodesics [6], [7]: the response of a detector depends on the 
frequency of a GW. 

The key point for an interferometric detector of GW is to measure the variation of 
the phase difference between two beams of light. The variation of the phase difference 
comes from the variation of proper lengths of two arms induced by a GW. On one hand, 
GW signals are submerged in the sea of noises and one has to find ways to extract the 
GW signals from noises. On the other, a detector needs the calibration to determine the 
sensitivity of the detector. Since both a GW and the mirror’s swings will result in the 
variations of the proper length of each arm, a naive idea is that the effect of GW on a 
detector can be mimicked by the vibrations of the mirrors’ position and many detector 
investigations have been made in this way [9], [10], [11], [13], [14]. The other purpose 
of the paper is to show that there is essential difference for the motion of photons in 
the background of GW and in a vibrating arm. The response of a detector depends on 
the frequency of a GW. For some specific frequency, the detector has no response. In 
contrast, the response of a detector to the vibrations with different frequency is almost 
the same. 

The rest of the paper is organized as follows. In section 2, the motion of photons 
in both a flat space-time and in the background of GW is studied in the framework of 
geometrical optics. The trajectories of photons are derived from the Fermat principle 
and are shown to be the same as those along null geodesics in general relativity. The 
time difference for a round trip between in the background of GW and in a flat space- 
time is shown to be dependent on the frequency of GW. In section 3, the propagation 
of the electromagnetic wave is investigated in wave optics. The deformed Maxwell 
equations and thus the wave equations for electric field are presented in the background 
of GW in a flat space-time approximation. The time difference is also calculated in 
the framework of wave optics. The results are the same as those in geometrical optics. 
Section 4 is devoted to the comparison of the differences between motion of photons in 
the background of GW and in a vibrating arm. The conclusions and remarks are given 
in section 5. Throughout the paper, the c = 1 units are used in the derivation and thus 
x =t. 


2. Motion of photons in the background of gravitational wave 


We are interested in the Minkowski spacetime and the spacetime with a plus-polarized, 
linear, normally incident, plane GW. In a time-orthogonal coordinate system, the 
metrics for the both cases can be written as 


ds? = —dt? + gijda‘de’, i,j =1,2,3. (1) 


The distance between two near points in a simultaneous hypersurface is measured by 
[16] 
dl? = gijdz’dat. (2) 


We consider a photon travels from point A to point B in space. There are infinitely 
possible classical paths. The length for each path is given by 


a dat dxi\'/ 


where dz'/dt is the coordinate velocity of the photon. 
In geometrical optics, the path of a photon is determined by the Fermat principle, 
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Here, only the isochronal variation is considered, so that ôt = 0. In the last step, the 
total-derivative term is absent because the end point of the integral is fixed. From the 
above expression, it is easy to see that 
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For convenience, we take l as the path parameter. Then, Eq.(5) becomes 
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In terms of the Christoffel symbols, it can be recast into the spatial part of 4-dimensional 
geodesic equation: 
da pde dr" 
TT 
It should be noted that Eq. (8) is different, in general, from the geodesic equation in 
3-dimensinal space, 


= 0, u,v = 0,1,2,3. (8) 


d's... de dz 
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The difference comes from the time-dependence of g;; in (8) or, in other words, from 


the time-dependent ‘refractive index’. 
Since the light travels along a null curve and since the path length should take 
extremal value, the traveling time of a photon from A to B 
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should also take the extremal value for fixed path. Namely, 
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Here, only the time variation is considered (6x = 0). In the second step, the total 
derivative is zero, too. Hence, 
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in terms of the Christoffel symbols, Eq. (12) reduces to the temporal part of four- 
dimensional geodesic equation: 


Ct. o dr” dr” 
dl? ” dl dl 


=0. (14) 


Eqs. (8), (14), (13) show that in the framework of the geometrical optics a photon 
travels along a null geodesic of a 4-dimensional space-time in both cases of a flat space- 
time and the background of GW, as expected in general relativity and the arc length 
l along a path serves as the affine parameter of the null geodesic. The trajectories of 
photons are the straight line in the 3 space just as the absence of a GW. The distance 
that a photon travels is the arc-length of the geodesic, L = f dl. Its value and thus the 
traveling time will be influenced by the time-dependent metric g,, in the background 
of GW. 

Now, we calculate the traveling time of a photon travels along x axis from mirror 
A at position x; = 0 to mirror B at zə = L. (From now on, L is the coordinate 
length of the arm.) When there is no GW and matter, the metric is of the form 
Juv = diag(—1,1,1,1). When a monochromatic GW with plus-polarization incidents 
along the z axis, the metric has the form 


guv = diag(—1, 1 + hi, 1 — hu, 1), hi = Rene C7], (15) 


where e11 is the amplitude of the GW, R means that the real part of the expression is 
taken, ¢ is the initial phase of the GW when the photon leaves mirror A. 
It is well known that if there is no GW, the distance between two mirrors in one 


ps / “ae (16) 


It is also the traveling time of a photon from one mirror to the other. In other words, 
Peak, 
In the background of GW, the traveling time of a photon from one mirror to the 


arm of the interferometer is 


other becomes T” and 


L L 1 
T' = f (g1) "dz = | (1 + sin) dz. (17) 
0 0 


In the last step of (17), the Taylor expansion has been used and the GW wave is 
supposed to be an extremely small perturbation around the flat space. Obviously, 


1 L 
0 


As expected, a GW background will change the length of the interferometer’s arm. In 
other words, a photon will ‘feel’ the distance difference AL = AT’ = (1/2) i hii (t)dax 
due to the presence of a GW. One should remember that hj, is a function of time t. 
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Now we consider a photon moves a round trip from mirror A to mirror B and then 
back to mirror A. For A — B, the distance that a photon travels is 


Th +B =T+ AT) ’ (19) 


and for B + A, the distance that a photon travels is 


0 1 tot+2T 
L to+T 
where to is the time when the photon leaves A. In the first line of the equation, dx 
is always negative. In the second line, the higher-order terms have been neglected. It 
should be noted that AT, # ATga in principle. 
For a round trip, we have 
1 


to+2T 
ar =z f haa GPS 
2 to 
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It shows that at specific values of the frequency of GW, f = 3 (with n = 1,2,3,...), 
AT’ = 0 no matter what value the initial phase ¢ takes. It means that for the GW with 
the frequency 5, there will be no measure difference between the presence and absence 
of a GW for traveling photons. The result has been reached in the literature [17], [8]. 
Define the dimensionless frequency of GW: a = 2Lf (> 0). Then, the traveling time 


of a photon for a round trip from A to B and back to A is 


AT’ = cubs, Jie") (i= a , (22) 


Ora 


In particular, when tọ = 0, 


AT = cue (sin (2ar + ¢) — sin (¢)] . (23) 


3. Propagation of electromagnetic wave in the background of gravitational 
wave 


In wave optics, light is treated as electromagnetic wave, which satisfies the Maxwell 
equations. 
The Maxwell equations in vacuum of a curved space-time read 


we pep) 
Ft = T= (VaR), = 0, (24) 


Favà F Ryu F Fusu =0 , (25) 


where Fw = —F,, = 0,A, — O,A, is the electromagnetic field tensor, A, is the 
electromagnetic potential. The electric field is 


E; = — F; U” , (26) 
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and the magnetic field is 
1 o 
B; = -3w Salve ’ (27) 
where e is the Levi-Civita tensor and U” is the 4-velocity of an observer, satisfying 


Juv U"U”" = —1. For static observers, U” = (1,0,0,0). Then, the electric field Æ and 
magnetic field B observed by the static observers are 


E; = — F; U” = — Fio = oÁ; = o; Ao i (28) 
1 NOT TI 1 Ao 1 ad „bo 
B; = T 5 Sim U" Py, = 9° Pio = 3I g  CiaB (OA, = OoAy) ` (29) 


In the background of GW (15), the explicit forms of three components of the 
magnetic filed B; are, 
Bı = (1 + h4) (32A; — 0342) , 
Bo = (1 — hy) (03A1 — 143) , (30) 
B; = (0,A2 — 02A)) . 


Therefore, 
Fo; = Fi 
Fiz = B3 , 
31 
Fə = (1 — hy) Bi, 3 
Fzı = (1 + h4) Bo, 
and, 
Fo = — (1 — h) Ei, 
PM = — (1+ h4) Bo , (32) 
Pi — E3 , 
F? = Bs , 
F” = By , (33) 
Ft — Bə 
In the background of GW (15), the Maxwell equations (24) read 
F o+ F”; =0. (34) 
Namely, 
F”, =F”; =0, (35) 
FY, +F*,=0, (36) 
or 


(hy — 1) Fit = (1 + hs) E22 = E33 =0 , 
((1 — h4+)E1) o + B32 — Bos =0, 
(1+ hy) £2) o — B3i+Bi3=0, 


Ezo + B21 — B12 =0. 
The Maxwell equations (25) 


Foi + Fjona + Fijo =0, 


Fije + Frij + Fini =0 , 


reduce to 
Ey 2 — E21 + B30 = 0, 


Eis — Eza — ((1 + h4) B2)o =0, 
E23 — E32 + ((1 — h4) Bı)o=0 , 
B33 + (1 + h+) B22 + (1 — h4) Bi =0. 


The above equations can be rewritten into a compact form 


where E and B are electric field and magnetic field vector in column matrix form and 
h is a square matrix: 


hn O 0 
h=| 0 -h 0], (38) 
0 0 0 


V is the gradient operator in 3-dimensional flat space. These are the deformed Maxwell 
equations in the background of GW in a flat space-time approximation. 

We try to setup the electromagnetic wave equations from the deformed Maxwell 
equations. First, making a time derivative at both sides of the second equation (37) and 
using the third equation, we get 


O(h- B) 
ðt ` 


The ratio of the frequency of GW and that of photons is very small and can be used as 
a perturbation parameter. Then, 


a vA E 


(39) 
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In components, the deformed electromagnetic wave equations are of the form 


OE 
(1 — hy) aa — VE, = —hy (Eru — Eza) — hy (F133 — Esa), 

0? By J 
(1+ hy) a2 VE = —hy (E112 — E222) — hy (E323 — E233) , (42) 
0° E 


OP V7? E3 = hy (E322 — E311) - 
The similar wave equations for B can be obtained in the same way. 

When the electromagnetic wave propagates along x axis, LE; = 0. Then, the first 
equation of (42) is the identity 0 = 0 in the leading order. The latter two equations 
read AE 
a = V? E =T h} E222 = hy (E323 = Ez 33) , 

(43) 
— V?E3 = hy (E322 — E311) - 


(1+ hy) 
0? E 
Ot? 
Since the right hand side of the first equation of (43) is 


= hy E222 + h4 (E222 + Ez 33) 
= h} E2 22 + hy (V7 E — E211) ; (44) 


the first equation of (43) reduces to 


0? Ey 


(1+ hy) PA 


(1+ hy) V? Eo = hi (E222 — Fan) . (45) 


The leading order of the perturbation of the electromagnetic wave’s equation is 


O° Ey 


3L — WES = h} (E222 — Ez 11) . (46) 


For the electromagnetic wave propagating along the x axis, E222 = 0 , E2833 = 0, 
E3 22 = 0 5 E3 33 = 0. Hence, 


or 

T = (1 — h,) Ez 11 = 0 . (47) 
The similar deduction gives rise to 

OE 

or — (1 a hy) E344 = 0 . (48) 


(47) and (48) are just the deformed form of the wave equation for electric field. The 
equations (47) and (48) gives the coordinate speed of the electromagnetic wave in the 


background of GW, 
h 
p= V1-hy 0al (49) 
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Similarly, one can obtain that 


w= Viti 14, (50) 
These mean that the photon’s coordinate speeds vary with the perturbation of space- 
time, which is caused by GW. The similar result has been obtained in other ways (say, 
in [8]). 
Suppose that the electromagnetic wave spend T’ and T to propagate from A to B 
and then back to A when there is and there is no GW, respectively. The time difference 
AT’ due to GW satisfies 


to+2T+AT’ to+2T+AT’ 1 
2L = | vdt = / (1 — lenene) dt. (51) 
t t 


0 0 


Again, set f = ṣ and tọ = 0. Eq.(51) results in 


enk (sin (ae 2b) 6) — sin 6) . (52) 


Its leading terms are again given by Eq. (23) because AT” is much smaller than L. 


Ar’ = 


Yar 


Obviously, when the dimensionless frequency a = 1, 2,3... , AT’ = 0, which means that 
the detector has zero response. 


4. Gravitational wave vs swing of end mirror 


In the above, only the ideal cases are discussed. Practically, an interferometer will 
undergo various vibrations. Whether it is possible to distinguish the GW from the 
swing of the mirrors is an important problem in the detection of GW. In the following, 
we shall show the difference between the effect of GW and the swing of a mirror. 

For simplicity, we consider the case that the mirror A is fixed and the mirror B 
undergoes a harmonic oscillation R[Ae~“?"/*+®)], where f is now the frequency of the 
vibration and ¢ is the initial phase of the B’s oscillation when the light leaves the mirror 
A. Let 7 be the traveling or propagating time of light for one-way trip when the mirror 
B undergoes an oscillation. Then, 


27 — 2L = R [2 (Ae *Prsete)y) (53) 
Denote the difference of the traveling or propagating time for the light’s one-way trip 


with and without end mirror’s swing by AT. We let a photon leaves mirror A at t = 0 
and consider the photon with initial phase ¢. Then, in the leading approximation, 


AT = Rl Ae CTIA) (54) 
AT = 0 requires that 
Rje tOr L+) =0, (55) 
1 
nfito=(S+n)q, (56) 
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Figure 1. Time difference of photon for a round trip from A to B and back to A 
as a function of frequency of GW. Parameters are chosen as to = 0, c = 3 x 10°m/s, 
e11 = 1072 and L = 4000m. The time difference also depends on the initial phase of 
GW. The different color curves represent the different initial phases. 


Obviously, the initial phase of the mirror’s swing will influence the difference of the 
traveling time. Since a continuous laser is used, the initial phases are uniformly 
distributed, equivalently. In other words, no matter what frequency the mirror B swings 
at, there are always some photons who can ‘feel’ the swing. In contrast, the photons 
cannot ‘feel’ the GW with frequency f = zF- 

Figure. 1 presents the plot of AT” due to a GW, given by Eq. (23). (The parameters 
are chosen as €;; = 107% and L = 4000m.) The different color denotes the different 
initial phase ø of GW. In the figure, the 8 curves correspond to @ = 0, 7, 5, an T, on an 
and a respectively. From the figure, we can see that the time difference of photon 
motions always vanishes when the frequency of a GW is 5. It means that a photon 
spends the same time to travel for a round trip in the background of specific GW with 
frequency >z as in a flat space-time. From Figure.1, we can also see that the initial 


phase of GW affects the readout in the measurement of GW. The envelop has the form 
of sinc z = sine 

Figure.2 shows the dependence of time difference of photon motions on the initial 
phase of GW for specific frequencies. It should be noticed that the curve corresponding 
with f = 4 (i.e. a = 1) coincides with the abscissa axis exactly. It means that the 
photon cannot ‘feel’ the GW with f = + no matter what the initial phase is. 


In the case of the mirror’s swing, the time difference of photon motions is 
AT = Acos (2rf L+ ¢). (58) 


We plot the frequency dependence of the time difference of photon motions caused by 
the mirror’s swing in Figure.3. Form the picture, we can see that the responses of 
the interferometer to different frequencies are almost the same. In this case, the time 
difference of photon motions also depends on the initial phase when this photon leaves 
the mirror A. The initial phase dependence is shown in Figure.4. 
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Figure 2. The initial phase dependence of the measurement for GW. The frequency 


of GW takes as f = 57, and a = 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6. 
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Figure 3. Frequency dependence of the time difference of photon motions caused 
by the mirror’s vibration. Parameters are taken as L = 4000m, c = 3 x 108m/s, A = 


4x 107!"m, again. For specific values of the initial phase, ¢ 
are depicted in different colors. 
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Figure 4. Initial phase dependence of the detector’s response. Parameters are taken 
as L = 4000m, c = 3 x 10°m/s, A = 4 x 107!"m. The different vibration frequencies 


are depicted in different colors. 
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5. Conclusions and remarks 


In this paper, the photon motions in the background of GW are re-analyzed in both 
geometrical and wave optics. In geometrical optics, the path of photons in 3-dimensional 
space is obtained by using the Fermat principle. The geodesic equations are re-produced 
and the photon paths are just the null geodesics in general relativity. The proper 
distance between the two mirrors and thus the traveling time that photons move in- 
between change with GW. In wave optics, the deformed Maxwell equations and thus 
the wave equations for electric fields are presented. The wave equations show that the 
coordinate speed v(t) of an electromagnetic wave moving between the two mirrors is 
modified by GW. The two methods give the exactly same equation (23) for the time 
difference influenced by GW. 

Both GW and the mirror swings result in the variations of the proper distance 
between the two mirrors. However, the characteristics of the variations are different. 
For a gravitational wave, the amplitude of AT” varies with the frequency of the GW, 
which is given by the envelop in Figure. 1. The property of the envelope has been 
reported in literature (see, for example, [6], [7], [8] etc.). For the mirror swing, however, 
the amplitude of AT is independent of the frequency of swings. The difference can be 
easily seen from the figures. The reason of the difference is as follows. A swinging mirror 
only affects the reflecting point of light, but does not affect the propagation of light in 
the arm. On the contrary, a GW affects the propagation of light in the arm everywhere. 
Therefore, the swinging mirror cannot mimic a gravitational wave very well, which is 
used in the previous studies on the instruments and calibration [8], [9], [10], [11], [12], 
[13], [14]. 

The above discussions are just for the motion of photons in a single arm. All 
the discussions are also valid for the motion of photons in the other arm. The only 
difference is that the phase variation in the two arms have different polarity due to 
the quadruple nature of the GW. The total effect should be the sum for two arms. If 
the time differences in both arms vanish for some specific frequency of GWs, then the 
detector has no response to these GWs. For the swinging mirrors, the detector may 
provide zero-response only when the two arms vibrates in synchronization. 

Finally, the time difference in a Michelson interferometer is discussed in the present 
paper. The detectors, like LIGO, Virgo, are all Michelson-Pabry-Perot interferometers. 
The analysis for a Michelson-Pabry-Pérot interferometer will be discussed elsewhere. 
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